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Abstract: Two independent Poisson streams of jobs flow into a single-server service system having 
a limited common buffer that can hold at most one job. If a type-i job (i = 1, 2) finds the server 
busy, it is blocked and routed to a separate type-i retrial (orbit) queue that attempts to re-dispatch 
its jobs at its specific Poisson rate. This creates a system with three dependent queues. Such a 
queueing system serves as a model for two competing job streams in a carrier sensing multiple access 
system. We study the queueing system using multi-dimensional probability generating functions, 
and derive its necessary and sufficient stability conditions while solving a boundary value problem. 
Various performance measures are calculated and numerical results are presented. 
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Files d'attente avec reemissions des clients et un serveur unique 



Resume : Un serveur unique 'exponentiel', sans file d'attente, est alimente par deux flux poissonniens 
et independants de clients. Un client du flux i (i = 1,2) trouvant le serveur unique inoccupe est servi 
et quitte le systeme une fois servi; s'il trouve le serveur unique occupe alors il rejoint une file d'attente, 
appelee orbite i, de capacite illimitee et equipee d'un serveur 'exponentiel'. Un client quittant l'orbite i 
est route vers le serveur unique; si celui-ci est inoccupe alors il quitte le systeme une fois servi; s'il est 
occupe alors il est de nouveau route vers l'orbite i et le processus decrit ci-dessus se repete. Le tout cree 
ainsi un systeme compose de trois file d'attente fortement correlees. Ce systeme de files d'attente sert 
a modeliser des protocoles de communication a contention. Nous etablissons une equation satisfaite par 
la fonction generatrice jointe du nombre de clients dans les trois files d'attente, que nous resolvons par 
reduction a un probleme aux limites de Riemann-Hilbert. Nous determinons la condition de stabilite du 
systeme et concluons le rapport en presentant des resultats numeriques pour les principales mesures de 
performance. 

Mots-cles : Files d'attente reemissions des clients, Probleme aux limites de Riemann-Hilbert, Protocoles 
de communication a contention 
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1 Introduction 



Queues with blocking and with retrials have been studied extensively in the literature (see e.g. [I]-[Z], 
[E]-[I1]> [IS]! HZ], [IS] and references therein). In this paper we investigate a single-server system with 
two independent exogenous Poisson streams flowing into a common buffer that can hold at most one 
job. If a type-i job finds the server busy, it is routed to a separate retrial (orbit) queue from which jobs 
are re-transmitted at a Poisson rate. Such a queueing system serves as a model for two competing job 
streams in a carrier sensing multiple access system, where the jobs - after a failed attempt to network 
access - wait in an orbit queue |231 124) . The two types of customers can be interpreted as customers 
with different priority requirements. An important feature of the retrial system under consideration is 
a constant retrial rate. The constant retrial rate helps to stabilize the multiple access system [B]. The 
retrial queueing systems with a constant retrial rate and a single type of jobs has been considered in 
H3]-[I1], HZ]- We formulate this system as a three-dimensional Markovian queueing network, 
and derive its necessary and sufficient stability conditions. Recently these stability conditions have been 
shown by simulations to hold even in a more general system with generally distributed service times [7]. 

The structure of the paper is as follows: After the Introduction we present the model in Section 2. 
Balance equations and generating functions are derived in Section 3, while necessary stability conditions 
are obtained in Section 4. Using the technique developed by Fayolle and Iasnogoroski [TH], in Section 5 we 
show that these generating functions are obtained, in closed-form, via the solution of a Riemann-Hilbert 
boundary value problem. This approach allows us to show that the necessary stability conditions found 
in Section 4 are also sufficient. Performance measures are calculated in Section 6, and numerical results 
are presented in Section 7. In particular, our numerical results demonstrate that the proposed multiple 
access system with two types of jobs and constant retrial rates provides incentives for the users to respect 
the contracts. 



2 Model 



Two independent Poisson streams of jobs, Si and S%, flow into a single-server service system. The service 
system can hold at most one job. The arrival rate of stream Si is Aj, i = 1, 2, with A := Ai + A2. The 
required service time of each job is independent of its type and is exponentially distributed with mean 
l/fi. If an arriving type-i job finds the (main) server busy, it is routed to a dedicated retrial (orbit) 
queue that operates as an -/M/l/oo queue. That is, blocked jobs of type i form a type-i single-server 
orbit queue that attempts to retransmit jobs (if any) to the main service system at a Poisson rate of Hi, 
i = 1,2. Thus, the overall system is comprised of three queues as depicted in Figure Q] 
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H l 






M-2 

Figure 1: Retrial system with two orbit queues. 

3 Balance equations and generating functions 

Consider the system in steady state. Let L denote the number of jobs in the main queue. L assumes the 
values of or 1. Let Qi be the number of jobs in orbit-queue i, i = 1, 2. The transition-rate diagram of 
the system is depicted in Figure [5J The numbers or 1 appearing next to each node indicate whether 
L = or L = 1, respectively. 

Define the set of stationary probabilities {P m „(fc)} as follows: 

Pnm(k) = P(Qi = m,Q 2 = n,L = k), m, n = 0, 1, 2, ... k = 0, 1. 
Define the marginal probabilities 

oo 

P m .(k) = J2 p mn(k)=P(Qi = m,L = k), m = 0,1,2,... fc = 0,l, 

and 



n=0 



P»n{k)= J2P mn (k) = P(Q 2 = n,L = k), n = 0,1,2,... k = 0,l. 



m=0 



Let us write the balance equations. If Qi = 0, we have 



(a) for Q\ — and k = 0, 

(b) for Qi = m > 1 and k = 0, 

(c) for Qi = and k = 1, 



APqo(O) =/iPoo(l), 
(A + /ii)P ro0 (0) =/xP m0 (l), 



(A + Ai)Poo(l) = APo (0) + miPio(O) + /i 2 P i(0), 



(1) 
(2) 

(3) 
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(d) for Qi = m > 1 and k = 1, 

(A + ft)Pmo(l) = AP mO (0) + /ilP m+ l, (0) + M2^ml(0) + A 1 P m _ 1 , (l). (4) 

If Qi = n, n > 1, we have 

(e) for Qi — and k = 0, 

(A + M2)Pon(0) =/xPon(l), (5) 

(A + //i+M 2 )P ro „(0) =/iP mn (l), (6) 



(f) for Qi = m > 1 and /c = 0, 



(g) for Qi = and k = 1, 

(A + /tx)P „(l) = AP „(0) + MiPm(O) + /i 2 P o ,„+i(0) + A 2 P ,„_i(l), (7) 

(h) for Qi = m > 1 and k = 1, 

(A + /i)P„ m (l) = AP„ m (0) + AtlP m+ l, n (0) + // 2 P m ,„+l(0) 

+AiP m _i,„(l) + A 2 P m , n _i(l). (8) 
Let us define the following Probability Generating Functions (PGFs): 

oo 

G { n ] (X)=Y / P ?nn(k)x m , k = 0,l, 71 > 0. 

rn=0 

Then, for n = and k = 0, multiplying each equation from (JTJ) and (J2|) by x m , respectively, and summing 
over m results in 

oo oo oo 



or 



(A + ^)G^(x) - miPoo(O) = (iG£\x). (9) 
Similarly, for n = and k = 1, using equations ([3]) and leads to 

oo oo 

(A + f i)G { a 1] = XG ( 0} +LiiJ2 Pm+i,o(0)x m + v 2 G ( ?\x) + X 1 Pn-i,o(l)z m . 

m— m— 1 

That is, 

(A + fi)G^\x) = XG^(x) + ^(G$\x) - Pqo(O)) + ^G^(x) + X lX G^(x). 
Multiplying by x and arranging terms, we obtain 

-(Xx + [i 1 )G i °\x) + (Ai(l -x) + A 2 +m)xG^(x) - M2.tG< 0) (x) = -/iiP 0O (0). (10) 
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Using equations §5§ and for n > 1 and k = results in 

(A + fx 2 )G^(x) + Mi(Gi o) 0r) - P o „(0)) = nG^(x), 

or 

(A + mi +M 2 )Gi 0) (a ; ) - MGi 1 '^) = Mi^bn(O). (11) 
Similarly, for n > 1 and fc = 1, equations Q and © lead to 

(A + = AG[°)(x) + ^(Gf (s) - P 0n (0)) + ^(x) 

+A 1 .TGW(x)+A 2 Gi 1 2 1 ( a; ), 

or 

-(Ax + /ii)G<?>(aO + (Ai(l - i) + A 2 + ^(i) - m xG { ^ +1 {x) 

-\ 2 xG i n 1 l 1 (x) = -/iiP o „(0). (12) 
Define now the two-dimensional PGFs 

oo oo oo 

J ffW(x, 2/ ) = ^^p m „(fc)x"Y i = E G « fc) (- T )y"' fc = ' L ( 13 ) 

n— ??i—0 n— 

Using equations ^ and multiplying respectively by y n and summing over n, we obtain 

(A + W )fr<°> (.*, j,) + M2 (7J(°) (x, j,) - G< 0) (a;)) - fiH^ (x, y) = ^H^ (0, y). (14) 
Similarly, using equations pop and (fl"2"]l . we obtain 

-(Ax + Hi)H (0) (x, y) + (Ai(l - x) + A 2 + fi)xH^(x, y) 

_^ (H (o) (a , )2/) _ G (0) (aj)) _ Aa^W^.y) - -/H# (0) (0, y). (15) 

Noting that Gq ^) = H^\x, 0) and denoting a := A + /ii + fx 2 , we can rewrite equations (fT4"| and (fT3|) 

as 

aHW(x,y)-»HW(x,y) = fx 2 H^(x,0) + ^H^^y), (16) 
(Axy + fi lV + fi 2 x)H^(x,y) - (X 1 (l-x) + X 2 (l-y)+^)xyH^{x,y) 

= fjL 2 xH^(x,0)+fM iy H^(0,y), (17) 

or, equivalently, in a matrix form 

C(x,y)H(x,y)=g(x,y), (18) 



where 

a —fx 
Xxy + fjL t y + [i 2 x -(Ai(l - x) + A 2 (l - y) + fi)xy 



C(x,y) = 
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U(x,y) 



H^(x,y) 



fi 2 xH^{x,0) + n iy HW(p,y) 



Now, if we calculate H^°\x, 0) and H^(0, y), the two-dimensional PGF H(x, y) is immediately obtained 
from equation (fTS|) . 



4 Necessary stability conditions 

Proposition 4.1 



ffW(l,l) = P(L = l) = - (19) 



and 



ff (0) (0,l) = P(Q 1 = 0,L = 0) = l-~ + — ) (20) 

/i \ /ii/ 

i/ (0) (l,0) = P(Q 2 = 0,L = 0) = l-~ + — V (21) 

The identities (fT?l) - (|2T|) show that conditions (i) A//i < 1 and (ii) (A/yu) (1 + Xi/fii) < 1 for i = 1,2, are 
necessary for the existence of a steady-state. Note that (i) is a consequence of (ii) so that in the following 
we will not consider condition (i) but only conditions (ii). 

Proof of Proposition l4TTt For each m = 0, 1, 2, ... we consider a vertical "cut" (see Figure [2]) between 
the column representing the states {Q\ = m,L = 1} and the column representing the states {Q\ = 
m + 1,L = 0}. According to the local balance equation approach [TT], we can write the balance of rates 
between the states from the left of the cut and the states from the right of the cut. Namely, we have 

AiP m .(l) = /XiP m+ i.(0), 771 = 0,1,2,.... (22) 
Summing (|22|) over all m results in 

Aiff«(l, 1) = - 1) - Po.(0)). (23) 

Clearly, P(L = k) = £^ =0 P m .(k) = i?«(l, 1), fc = 0, 1. 
From (|23[) we readily get 

l_p o .(0) = ^±a f r(i)(i,i). (24) 

Ml 

Since P o .(0) = H^°\0, 1), we can write (JMJ) as 



l-£T<°)(0,l) = ^±^ff«(l,l), (25) 
Mi 
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and, by symmetry, 

l-tf^(l,0) = ^±^ff«(l,l). (26) 
Substituting (|2"!j)) and (|26|) in equation (fTB")) , with x = y = 1, yields 

= P(Z = 1) = -. 

Now, from (|25| and (|26| . respectively, we obtain 

ff<°>(0,l) = P(Qi = 0,L = 0) = l-- f^l±i^ N ) (27) 



M V Mi 

and 



(1,0) = P(Q 2 = 0, X = 0) = 1 - - ( ) , (28) 



which completes the proof. I 

The next result shows that the system cannot be stable if either (A//i)(l + Ai/ fix) = 1 or (A//i)(l + 
A2/M2) = 1. 

Proposition 4.2 If either (X/fi)(l + = 1 or (A//i)(l + X 2 /fi 2 ) = 1 then P OT ,„(0) = P m ,„(l) = 

/or a?/ m, n = 0, 1, . . . or, equivalently, both queues Q\ and Q2 are unbounded with probability one. 

Proof. Assume, for instance, that (A//x)(l + A2//12) = 1 so that H^°'(1,0) — from (|2T|) . Since 
iP°)(l,0) = Em>o p ™,o(0) (see P^jl). the condition iP°)(l,0) = implies that 

P m ,o(0)=0 for m = 0,1,..., (29) 

so that from CD"© 

Pm,o(l)=0 form = 0,1,... . (30) 

We now use an induction argument to prove that 

P ro ,„(0) = form, re = 0,1,.... (31) 

We have already shown in (f2"9"|) that (f5Tj) is true for n = 0. Assume that (|3T|) is true for n = 0, 1, . . . , k 
and let us show that it is still true for re = k + 1 . 

From ([6]) and the induction hypothesis we get that Pm,fc(0) = P m ,fc(l) = for m — 1,2,.. .. The latter 
equality implies, using ©, that P m) fc + i(0) = 0. This shows that (|3"Tj) holds for m = 0, 1, ... and n = k+ 1, 
and completes the induction argument, proving that (|31D is true. 

We have therefore proved that P m ,„(0) = for all to, n = 0, 1, . . .. Let us prove that P mi „(l) = for all 

m,n = 0, 1, The latter is true for m,n = 1,2,... thanks to ©. It is also true for n = 0, m = 0, 1, . . . 

from pop . It remains to investigate the case where to = and n = 0, 1, By ([5]) and (|3"Tj) we get that 

Po,n(l) = for re = 1, 2, . . ., whereas we have already noticed that Pq ! o(1) = 0. 
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In summary, P m , n (0) = P m .n(l) = for all m, n = 0, 1, . . ., so that 

P(Ql = m, Q 2 = n) = Pm,n(0) + Pm,n(l) = 

for allm,n = 0,l,..., which completes the proof. I 
We conclude from Propositions (|4.1[) and (|4.2p that conditions 

^Ui+M<l and (^]( l + h\ <l (32 ) 
M/ V Mi/ \M/ V M2/ 

are necessary for the system to be stable. We will show in Section [5] that these conditions are also 
sufficient, thereby implying that they are the stability conditions of the system. 



5 Derivation of (x, 0) and (0, y) 

Throughout we assume that the necessary stability conditions found in (p?2"T) hold. Our analysis below 
will formally show that these conditions are also sufficient for the stability of the system. Let us give an 
intuitive motivation for this result. In a stable system, A//i is the fraction of time the server in the main 
queue is busy. Thus, this is also the proportion of jobs sent to the orbit queues. Therefore, the maximal 
rates at which jobs flow into orbit queue 1 and into orbit queue 2 are (Ai + fj,\)X/n and (A2 + yU 2 )A/^, 
respectively. Each of these rates must be smaller than the corresponding maximal service rate, ^1 or // 2 , 
respectively. 

Lemma 5.1 Conditions iSty) imply that either a\i < or 01X2 < Ht l 2- 
Proof: Assume that aAi > /x/ii and a\ 2 > (J.fJ-2 

Multiplying the first inequality in (|32p by fifii and using the definition of A and a gives 

(Ai + A 2 )(Ai + /xi) < < aXi = (Ai + A 2 + fi± + ^2)^1 

which is true if and only if (a) X 2 [i\ < Ai//2- 

Multiplying now the second inequality in (|32[) by [i[i 2 gives 

(Ai + A 2 )(A 2 + H2) < MM2 < oiX 2 = (Ai + A 2 + fii + ^ 2 )A 2 

which is true if and only if (b) Ai/x 2 < A 2 /xi. 

Since inequalities (a) and (b) cannot be true simultaneously we conclude that either aAi < fifi± or 
aA 2 < fj,fj,2, which concludes the proof. ■ 

From equations (|16p - (|17p we obtain the two-dimensional functional equation 

R(x, y)H<® (x, y) = A(x, y)H^ (x, 0) + B(x, y)H<® (0, y), \x\ < 1, \y\ < 1, (33) 
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with 

R(x,y) 
A(x,y) 
B(x,y) 

For further use note that 



= Aia(l - x)xy + X 2 a(l - y)xy - £t//i(l ~ x )v - m(l - V)x (34) 
= ((i-y)(A 2 y-/i) + A x (l-x)y)n2X (35) 
= ((l-x)(Xix-fi) + \2(l-y)x)fi 1 y. (36) 



R{x,y) = — A(x,y) + - y)x + nm(x - y), (37) 

Oi 

R{x,y) = — B(x,y) + Xn(l - x)y + fj,/j, 2 (y - x). (38) 
Mi 

The kernel R(x,y) of the functional equation (j3"3"]) is the same as the kernel in |18l Eq. (1.3)] upon 
replacing Xi and fii in |18j by Xia and fj,ifi, respectively, for i = 1,2. 

In the following we set Xi = aXi and pn = for i = 1,2. In this notation, the kernel R(x, y) is expressed 
as 

R(x,y) = Ai(l - x)xy + A 2 (l - y)xy - /ti(l - - /i 2 (l - y)x. (39) 
Also define A = Ai + A 2 = aX. 

Assumption A: Without loss of generality thanks to Lemma 15.11 we will assume throughout that 
aAi < or, equivalently, that Ai < 

Once H(°\x,y) is known for all |x| < 1 and \y\ < 1 then H^\x,y) can be found from ([TBI . In the 
following we will therefore only focus on the calculation of H^°'(x,y) or, equivalently from (|3"3"]1 . on the 
calculation of H<® {x, 0) and H^(0,y) for all \x\ < 1 and \y\ < 1. 

We will show in Section IQ1 that if (°)(x, 0) is given by the solution of a Riemann-Hilbert problem on the 



circle centered at x = and with radius y /ti/Ai (see ([6^]) ). from which we will derive H^°\0,y) for all 
|y| < 1 (see 

The technique of reducing the solution of certain two-dimensional functional equations (equation (|33|) in 
our case) to the solution of a boundary value problem (typically Rieman-Hilbert or Dirichlet problem) - 
whose solution is known in closed- form - is due to Fayolle and Iasnogorodski [TS]. In [TB] (see also [2"0] 
that generalizes the work in |18| ) the unknown function is the generating function of a two-dimensional 
stationary Markov chain describing the joint queue-length in a two-queue system. Cohen and Boxma |15) 
extended the work in |18l [2Tj] to two-dimensional stationary Markov chains taking real values, typically 
representing the joint waiting time or the joint unfinished work in a variety of two-queue systems. Other 
related papers include Ell EE U3\ (non-exhaustive list) . 

5.1 Branching roots of R(x,y) 
For y fixed, R{x, y) vanishes at 



-6(g) ± y/c(y) 

x{y) = ~i (40) 

2X iy 
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where 



with 



We have 



with 



b(y) := X 2 y 2 - (£1 + h + \)y + A2 (41) 

c(y) := b-(y)b+(y) (42) 

6_(y) :=b(y)-2y^M, b + {y) :=b{y) + 2ysJ\^. (43) 

b-(y) = \ 2 (y - yi)(y - y 4 ), b+(y) = X 2 (y - y 2 )(y - 2/3) (44) 



" = — ^ ' m = (45) 



6 + - 4A 2 a 2 6 + v^? -4A 2 a 2 . . 

2/3 = — 2X — ' 2/4 = — 2X — (46) 



£1 = Ai + A2 + A + 2yAiAi, £ 2 = Ai +A2 +A-2yAiAi- (47) 

j/i, . • • , 2/4 are the branch points of x(y) (since c(yi) = for i = 1, . . . ,4). It is easily seen that (Hint: 
2/2 < 1 and 2/3 > 1, both from Assumption A)) 

< 2/1 < 2/2 < K 2/3 < 2/4- (48) 



Remark 5.1 The algebraic function x(y) has two algebraic branches, denoted by k(y) and k" (y), related 
via the relation k(y)k a (y) — A1/A1. When y £ (2/1,2/2) U (2/3,2/4) Ml/) an ^ k° (u) are complex conjugate 
numbers (since c(y) < for those values of y ), with k(yi) = k a (yi) for i = 1, ...,4. In particular, 

\k(y)\ = V k (y) ka (y) = \J A1/A1 /or y £ [2/1,2/2] U [2/3,2/4], thereby showing that for y £ [2/1,2/2] (resp. 

2/ G [2/3,2/4] ) k(y) arl ^ k a (y) lie on the circle centered in radius *J A1/A1 . 

When x is fixed similar results hold. We will denote by 



2/W = (49) 

2A 2 x 

the algebraic function solution of R(x, y) = for x fixed, where e(x) := Aiir 2 — (Ai + A2 + A)x + Ai and 
:= e_(a;)e+(a;), with 

e_(x) := e(x) - 2x^\~ip- 2l e + (x) := e(x) + 2x\j \ 2 (i 2 . 

We denote by Xi, i = 1, . . . , 4 the four branch points of y(x), namely, the zeros of d(x); they are obtained 
by interchanging indices 1 and 2 in (|4"5|) - (l4"Tl) . 
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We have 

e-(x) = \i(x — xi)(x — X4), e+(x) = \i(x — X2)(x — X3) (50) 

where 

< x\ < X2 < 1 < £3 < X4 (51) 

with X2 = 1 iff A2 = fi2 ■ 

The following results, found in |18l Lemmas 2.2, 2.3, 3.1], hold : 

Proposition 5.1 For y fixed, the equation R(x,y) = has one root x(y) = k(y) which is analytic in the 
whole complex plane cut C along [yi , 2/2] and [2/3 , 2/4] . Morecwej[i] 

(al) \k(y)\ < 1 if \y\ = 1. More precisely, \k(y)\ < 1 i/ |y| = 1 ure£/i y^\, and k(l) = min(l, /ti/Ai) = 1 
under Assumption A. 

fMJ |fc(2/)|< ,/f forallytC; 

V A i 

fcij w/ien ?/ sweeps twice [f/1,2/2]; k(y) describes a circle centered in with radius \J~^ , so that \k(y)\ 



A 



/or y g [2/1,2/2]- 



Similarly, for x fixed, the equation R(x, y) — has one root y{x) = h(x) which is analytic in<C — [xi, X2] 
[x3, 2:4], and 

(a2) \h(x)\ < 1 if\x\ = 1, x ^ 1, and h(l) = min(l, A2/A2) < 1. 
0>2) \h(x)\ < for all x 6 C; 

V ^2 

(c2) \h(x))\ = ^ ifxG[ Xl ,x 2 } 
Moreover, 

(dl) h(k(y)) = yfory£ [2/1,2/2] and (d2) k(h(x)) = x for x € [2:1,2:2]. 



(d2) h(\J '/ti/Ai) = 2/2 a^rf h(-y]h/\i) = Vi- 



(d3) k{\J A2/A2) = £2 and \/ ' ^2/^2) = x\. 
Last 

(e) \h(x)\ < 1 /or 1 < |x| < (VecafZ £/ia£ Ai < /ij. 

1 Apply Rouche's theorem to _R(x, j/) to get (al), and the "maximum modulus principal" to the analytic function k(y) in 
C — [?/i,S/2] — [2/3,2/4] to get (bl). (cl) follows from Remark l5.ll 
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5.2 A boundary value problem and its solution 

We are now in a position to set a boundary value problem that is satisfied by the unknown function 

In the following, C a = {z G C : \z\ < a} (a > 0) denotes the circle centered in of radius a, and 
C+ = {z G C : \z\ < a} denotes the interior of C a . 

We know that R(k(y),y) — by definition of k(y). On the other hand, H^°\x,y) is well-defined for 
all (x,y) = (k(y),y) with \y\ = 1, since (i) H^°'(x,y) is well-defined for \x\ < 1, \y\ < 1, (ii) k{y) is 
continuous for \y\ = 1 (from Proposition 1 5 . 1 1 we know that k(y) is analytic in C — [2/1,2/2] and we know 
that < 2/1 < 2/2 < 1 so that k(y) is continuous for \y\ = 1), (iii) \k(y)\ < 1 for \y\ = 1 (cf. Proposition 
IS-U fal)). Therefore, the l.h.s. of (|3"3")l must vanish for all pairs (k(y),y) such that \y \ — 1, which yields 

A(k(y), y )H^(k(y) 1 0)^~B(k(y),y)H^(0,y), V|y| = 1. (52) 

The r.h.s. of (|5"2")) is analytic for \y\ < 1 with y £ [2/1,2/2] and continuous for \y\ < 1, so that the r.h.s. of 
(|52p can be analytically continued up to the interval [2/1 , 2/2] ■ 

This gives 

A(k(y),y)H^(k(y),0) = -B(k(y),y)H^(0,y), \/y G [yi , 2/2] . (53) 

It is shown in Lemma [A. II that B(k(y), y) 7^ for 2/ £ [2/1, 2/2]- We may therefore divide both sides of (|5"3"|) 
by B(k(y),y) to get 

|||^^°)(fc(2/),0) = -i/(°)(0,2/), V2/G [2/1,2/2]. (54) 
Take 2/ G [2/1,2/2]: we know by Proposition 15. H -(cl) that k(y) = x G C / „ » so that h(k(y)) = y = 

v Ml / "l 

G [2/1,2/2] by Proposition 15. II (dl). We may therefore rewrite (|54")1 as 

||^| ff ™(,,0,^,/».(0,Mx)), V« C/ -. (55) 



It is shown in Lemma [A. 2 1 that h(x) is analytic for 1 < \x\ < yfii/Xi and continuous for 1 < \x\ < 



/ti/Ai; furthermore \h(x)\ < 1 for 1 < |x| < y/ti/Ai by Proposition 15. li fe). These two properties 



imply that, H t ^ ' > (0 1 h(x)) 1 the r.h.s. of ([55]) . is analytic for 1 < \x\ < yAi/Ai and continuous for 
l<\x\< yAi/^i . which in turn implies that, B^j g^^O), the l.h.s. of ([Ml), can be extended as 



a function that is analytic for 1 < \x\ < y jli/Xi and continuous for 1 < |x| < W/ii/Ai 



It is shown in Lemma I A. 31 that A(x, h(x)) has exactly one zero in (1, y /ti/Ai], of multiplicity one, given 

by 

-(A + /-ti - n)Xfii + y/((A + - m)^Mi) 2 + 4AAi(A + Mi)Wfi ft -^ 
X ° 2AAi(A + m ) ' (55j 



if x < y fii/Xi and if (A + /ii)xo/(Axo + /Xi) < y ^2/^2 and does not have any zero in (1, y /ti/Ai], 
otherwise. 
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Introduce 



A(x,h{x)) 



U(x) := ' v and ff(ar) := 0)(x - x ) r , (57) 

£>(x, n(x))(x — xq) 1 

where r G {0, 1} is defined by 

1, ifxo < J Ai/Ai and A±B_g° < J~fo/\ 2 

v Ax + ^i v (58) 

0, otherwise. 



By construction 



As noticed earlier the l.h.s. of (|5T))) is analytic for 1 < |ar| < y Ai/Ai and continuous for 1 < \x\ < y /zi/Ai 



(59) 



Since by construction t/(a;) does not vanish in (1, y /ii/Ai] we conclude from (|59p that the function H(x) 
that is initially analytic for |x| < 1 and continuous for \x\ < 1 can be extended as a function that is 



analytic for \x\ < y fli/Xi and continuous for |x| < y /Si/Ai. 
In summary, we have shown that the real part 

5ft (i U(x)H(x)j = 0, Vie C^—jj^, (60) 
where H(x) is analytic in C + , — and continuous in C + , — U C / — r*— , and where U(x) does not 

w J VaTaT _ _ \Wm VA1/A1' _ v ; 

vanish on C , — 7*-. This defines a Riemann-Hilbert boundary value problem on the circle C , — r=— , 
whose solution is given below. 

Define 

X := [arg £/(a;)] xeC (61) 

7T vAl/*l 

the so-called index of the Riemann-Hilbert problem, where [arga(z)] zS c denotes the variation of the 
argument of the function a{z) when z moves on a closed curved C in the positive direction, provided 
that a(z) ^ for z G C). 

The Riemann-Hilbert problem has x + 1 independent solutions [22} P- 104]. It is shown in Lemma TA.4I 
that, as expected, x = under conditions (|32p . thereby showing that the solution of the Riemann-Hilbert 
problem (|60p is unique under conditions (|32[) which will in turn imply that (|32p are sufficient stability 
conditions for the queueing system at hand. 

With x — the solution of the Riemann-Hilbert problem is 

H(°\x,0)=D(x-x )- r exp( -L / log(J(z)) dz ) , V|a?| < a/MiMi, (62) 

J|^| =a /a77a7 z-x J v 

where Z? is a constant (to be determined) and (with z the complex conjugate of z G C) 



J(,) = -^ 



Inria 



Files d'attente avec reemissions des clients 



15 



We are left with calculating the constant D in (JH2J) ■ Setting x = 1 in ([52)1 gives 

D = (l-x„r(l-i(l + ^)exp(--L/ MiMU") (63) 



by using the value of if(°) (1,0) found in (|2T|) . We may therefore rewrite (j62|) as 

ff <.» (l ,0 ) = (-i^»y( 1 -^( 1 + ^))e x pf-L/ ]«£_11^ (64 , 

\x-xqJ V Ma// V 2tti J| 2 | = y^7I7 (z-a:)(z-l) 



for all |a?| < y AiMi- 

We also need to calculate the other boundary function H^-°\0,y) for \y\ < 1. For |y| = 1, H(°\0,y) is 
given in ([5^)1 . For |y| < 1. H^(0,y) is obtained from ([5^1) and Cauchy's formula, which gives 



2tti J| t | =1 i - y 



Itfl < 1- (65) 



where 

^■■'-wm l Hmm -° l <66) 



does not vanish for all lil = 1, as shown in Lemma [AJ3] 



Introducing and (|55|) into (fT5|) uniquely determines the joint generating functions H^(x,y) and 
H( x >(x,y) for |x| < 1, |y| < 1 which shows, as announced, that conditions are also sufficient for the 
system to be stable. 



6 Performance measures 

Later on in this section we shall need the derivatives £H^(x,0)\ x=1 and j-# (o) (0>J/)ly=i- 
Differentiating (|64| w.r.t a; gives 

ax \ x — xo ) \ fJ. \ A f 2 

x eJ—l log^z))^-!) ; 

\ 27r W| 2 |=v^7l7 (z-x)(2-i) / 

I -r If log(J(z) \ 



2 



- H^, 0) f -=L_ + J- f ^^dz) . (67) 

\x-xo 2m J\ z \ = ^- 717 (« - x) 2 J 

Letting x — 1 in (j6~?| and using (j2~Tj) yields 

-f ff <.>( I ,„ )| „ I = ( 1 _^( 1+ M)('^ + 1 / *5fcM*Y (68) 
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The derivative -r-if( o )(0, y)\ y =i is obtained from ([52"]) . By Lemma [A. 51 we have 



dy 

where (x,0)\ x=1 and ff(°)(l,0) are given in (|68p and (|2ip . respectively. The limits in the above 

expression can be calculated by L'Hopital's rule. Lengthy but easy algebra gives 

A(k(y),y) _ (A 2 -/i + Aifc'(l))/i 2 



y^iB(k(y),y) (A 2 + (Ai - /x)fc'(l))/xi 
and 

y-n dyB(k(y),y) 

(-A 2 + (-Ax + ^)fc'(l) + (A 2 - ^)fc'(l) 2 + Axfc'(l) 3 + (/x — Ai — A 2 )fc"(l))/^ 2 

(A 2 + (Ai — ^i)fc'(l))/ii 

where 

A 2 - A2 



fc'(l) 

A'i - Ai 
and 

fc"(i) = 2-^ + '" 2 ~ + ^ 2 ^iA2 + ^i^ 2 + -^Ai 



(Ai-Ai) 3 

We are now in a position to calculate some important performance measures. 

By setting x = in equation (|64p . we immediately obtain the probability of empty system 

P(Q 1 = 0,Q a = 1 £ = 0)=(^YYl-^(l 



xexpf-L/ (70) 
V^W^/XT z{l-z) ) 

Next, we calculate the expected orbit queue lengths. For the first queue, we have 

OO / OO OO \ j 7 

E[Q 1 }= ^P„ m (0) + ^P m „(l) = _ff(°)(x,l)U =1 + — if (x,l)| x=1 . (71) 

m=l \n=0 n=0 / 

Thus, we need to calculate -^H {0 \x, l)| x =i and ^ijW(x, l)U=i. From ([33]) we have 

ff( ° )(s,, ' ) = W~y) H{ ° ){x ^ ] + W^) H{0){0 ' y) - (72) 

Using d3J])-(I3SJ) and setting y = 1 in ([72]). yields 

ff(o) (a . i) = g(0)( 0) + te^ g (o) (0i i). 
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Next, by differentiating the above relation with respect to x we get 

ax yaXix — pL^LiY aX\x ~ fifi\ ax 



(a\ix - im^Y 



Setting x = 1 in the above, yields 



Al " 2 d ffC'fx.OJU,,, (73) 



/Lt/Ui — aAi da; 

where ff(°)(0,l), ff(°)(l,0) and dH^(x,0)/dx\ x=1 are given in (|20p . (|2ip and (|68p . respectively. 

It remains to find dH^(x, l)/dx\ x =i- Differentiating §TE\i with respect to x and setting x = y = 1 gives 

±HV{x,l)\*=l = -^ (0) (x,l)U =1 -^-f^(o) (X;0)U=1 
ax /i ax [i dx 



5 -^ p ( (o -,, 1 )^>(0,l)-M 2 H<»)(l,0)) 



li -H<»'(x,0)|„ 1 , (74) 



///ii — aAi dx 



by using (|75|) . 

By combining ([TTjl. {75} and ([75]) we finally obtain 

/i/ii — aX\ dx 

where ff( )(0,l), ff(°)(l,0) and dH<-°\x,0)/dx\ x=1 are given in (|20p . (|2ip and (|68p . respectively. 
Similarly, the expected queue length for the second orbit is given by 

E[Q 2 ] = ^ (0) (l,»)l«=i + ^ (1) (l,i/)lv=i 



(75) 



(a + ^ )A fi - ^ 2 )^ (0) (1, 0) - Mi^ (0) (0, 1) 
«u 2 - ckAt ^ V 



(ju/x 2 - «A 2 ) 
_ M1 (A 2 + M 2 W (0) 
/i^ 2 - aA 2 ay 



(76) 



where dH^(0, y)/dy\ y =i is given in ([BUI) . 
Finally, we recall that (see (fTTJIO 



£[L] = P(L=1) = -. 

A 4 
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7 Numerical examples 

To obtain more insights into the performance of the system, let us consider numerical examples. First, 
we set = fj,2 = 2, /i = 4, Ai = 0.1 and vary A2 in the interval [0.2; 1.9]. In Figure [3] we plot the 
probability of an empty system P(Qi = 0, Q 2 = 0, L = 0) calculated by ([70)) as a function of A2. We also 
plot -ff(°)(l,0), see formula (f2~T)) . which corresponds, if Ai is small, to the probability of empty system 
with one type of jobs and a single orbit queue. Now if we change the value of Ai from 0.1 to 1.0, we 
observe that the value of P(Qi — 0, Q2 = 0, L = 0) deviates significantly from H^°\l, 0). 

Keeping \i\ — H2 — 2, fi = 4, in Figure [4] we plot the expected queue length of the second orbit E[Q2] 
calculated by ([75)) as a function of A2 for Ai = 0.01; 0.1; 1.0. We also plot the expected queue length of 
the orbit queue for the single orbit retrial system 0], which is given by 

prnl A^(A 2 + ^ + M2) 

|U(jU/X2 - A| - \2ll2) 

Again, as expected, when Ai goes to zero, £ , [Q 2 ] approaches the expected queue length of the orbit queue 
in the single orbit retrial system. 

Next, we investigate how the retrial rates affect the system performance. Let us fix Ai = A 2 = 1.2, \i = 4, 
/Lti = 2 and we vary fj,2 in the interval [2.0; 2.15]. With such parameter setting, the system is not too far 
from the stability boundary. We plot in Figure [5] the expected lengths of the orbit queues, i?[Qi] and 
E[Q2], as functions of ^2- We can see that if the jobs of type 2 retry at a bit faster rate than the jobs 
of type 1, they can gain significantly in terms of the waiting time. Specifically, an increase of less than 
10% of the retrial rate of jobs of type 2 helps them to reduce the expected orbit queue length by 50%. 
Clearly, if there is no cost for retrials, it is beneficial for the jobs to increase their retrial rate. However, 
there are good reasons to keep the control of the retrial rates in the hand of the system administrator 
and not to set them too high. As was just mentioned, the first reason is the possible cost for retrials. The 
second reason is the creation of incentives to respect the contract. To illustrate this point, we fix Ai = 1, 
/ti = /i2 = 2, fj, = 4, and vary A2 in the interval [0.2; 1.34]. In Figure [SI we plot the expected queue 
lengths of the orbit queues. We see that if the jobs of type 2 increase their input rate beyond their fair 
share, they will be severely penalized in terms of the expected delay, whereas the increase of the input 
rate of jobs of type 2 does not inflict any significant damage to the jobs of type 1. 
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A Appendix 

Lemma A.l Under conditions i3ty), (i) A(k(y) 1 y) 7^ and (ii) B(k(y),y) 7^ for 

y e [3/1,2/2]- 
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Equivalent]), (Hi) A(x, h(x)) 7^ and (iv) B(x, h(x)) ^ for x € C / . . 

Proof. From l[55]l and l[3"8]) we see that R(x,y) and B(x,y) vanish simultaneously if and only if 

(1 -x){\xx-n) + \ 2 {l-y)x = 
A(l - x)y + pt 2 {y - x) = 0. 

The second equation gives x = (A + fJ- 2 )y/ ((Ay + M2))- Plugging this value of 2 into the first equation 
yields (Hint: use A = Ai + A2) 

Pi(v) := (l-y)Qi(y) 
with Qi(y) := AA 2 (A + M2)y 2 + (A + M2 - m)Am22/ - MM 2 = 0- 

From lim 1) _j.± 00 Qi(y) = +00 and Qi(0) = — MM2 we conclude that the polynomial Qi(y) has two real 
roots, y_ < < y+ and that Qi(y) < for < y < y+. Since 

Q 1 (l)=(^)^(l + ^)-l)<0 > (77) 

where the latter inequality holds under conditions (|3"2"j) . we conclude that Qi(y) < for y € [0, 1], which 
in turn implies that Pi(y) < for y G [0, 1). The latter completes the proof of (ii) since [yi,y2] C [0, 1) 
(see (|15]>). 

The proof of (i) is the same as the proof of (ii) up to interchanging incides 1 and 2. 

Eqns (iii) and (iv) both follow from the fact that fcQyi, y 2 ]) = C ^ nj^j^ (cf- Proposition 15 . 11 (11)1 and the 
relation h(k(y)) = y for y <G [yi,y2] (cf. Proposition 15 . 11 fdl ) ) . ■ 



Lemma A. 2 Under condition A, h(x) is analytic for 1 < |x| < y M1/A1 and continuous for 1 < |x| < 
V^Ai/Ai 

Proof. We already know by Proposition lS.ll that h(x) is analytic for x € C—\x\, x 2 ] — [x^, X4] where x 2 < 



1 < X3 . It is therefore enough to show that y Mi / Ai < £3 or, equivalently from (|50)l that e+ y /ti/Ai I < 



0. Easy algebra shows that e + ( y M1/A1 ] = — y/ti/Ai ( ( vAi — + ( V A2 + VM2) ) < 0, which 

concludes the proof. 



Lemma A. 3 Assume that conditions h3ty) hold. Define 



_ -(A + mi - m)A/xi + y/((A + /ii - /i)A^i) 2 + 4AAi(A + mQmm| 

2AA!(A + mi) > 

If Xo < y^ A1/A1 and i/ (A + fj,i)xo/(Xxo + Ml) ^ y^ M2/A2 £/ien -A(x, /ias exactly one zero x = x$ in 
1, v/ M1/A1 and t/iis zero /ias multiplicity one. Otherwise A(x, h(x)) has no zero in ( 1, \J fii/Xi 
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Proof. From |55|) and (|57|) we see that R(x, y) and A(x, y) vanish simultaneously if and only if 

(1 - y){X 2 y - fi) + X^l - x)y = 
A(l - y)x + m(x - y) = 0. 

The second equation gives 

y= { -^>. (78) 

Ax + ll\ 

Plugging this value of y into the first equation yields 

(Xx + my 

with Q 2 (x) := AAi(A + Hi)x 2 + (A + [i\ - fi)XfiiX — /j.fj,f. 

The polynomial Q 2 {x) has exactly one positive zero given by xq. From the inequality 

Q 2 (l) = m(A + Mi) (- + — - < 



which holds from (p?2"1) . together with lim a; _j.± 00 Q2(z0 = +oo and Q2(0) < 0, we conclude that 1 < xq. 
This shows that 



If xq > y fli/Xi then A(x, h(x)) has no zero in (1, y jli/Xi]; 



Assume that xq < y (L\/\i. A(x, h(x)) as a unique zero in (1, y /ii/Ai], given by x = xq provided 



that (see (|75|) /i(a;o) = (A + fii)xo/(Xxo + Hi) < y A2/A2 since we know from Proposition 15 . 11 - (b2 s ) 



that the branch h(x) is such that \h(x)\ < y (J.2/X2 for all x £ C; if (A+/Ui)a;o/(Aa;o + /ii) > V A2/A2 



then A(s, /i(a;)) does not vanish in (1, y /ii/Ai]. 

We are left with proving that when ^4(x, ft.(x)) vanishes at x = xq then this zero has multiplicity one. 
From now on we assume that A(xq, h(xo)) = 0. 

From the definition of h(x) and (|57)) we get 

= R(x, h(x)) = — A(x, h(x)) + /dA(l - h(x))x + ui(x - h(x))]. 

Differentiating this equation w.r.t. x gives 

= — dA ( x > h ( x )) + fX [-Xh'(x)x + A(l - h(x)) + Mi(l - h'(x))]. (79) 
fi2 dx 

Assume that dA(x, h(x))/dx = Oat point a; = xq, namely, assume that A[x, h(x)) hasazero of multiplicity 
at least two at x = xq. From ([?9")l this implies 

-Xh'{x )x + A(l - h(x )) + /ii(l - ti(x Q ) = 
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that is 



h'(x ) = m ,, X+ J 11 , 2 (80) 
(Xx + HiY 



with h(x ) = (A + (1i)xq/(\x + /Lti) (see ([75]) . 

On the other hand, letting (x,y) = (x, h(x)) in (|35p yields 

h(x)) = ((1 - /i(x))(A 2 ft.(a:) - //) + Ai(l - x)h{x))pL 2 x. 

Differentiating ^4(x, /i(x) wrt x in (|8 If) and letting x = xq gives 
dA(x, h(x)) 

Tx u=xa = 

[-h'(x )(Xh(x ) - n) + A 2 (l - h(x ))h'(x Q ) - Xih(x ) + Ai(l - Xo)h'(x )]^ 2 Xo 

+—A(x ,h(x )) 
x 

= [h'(x )(-2X 2 h(x ) + A 2 + [i + Ai(l - x Q )) - Ai/i(x )]/i 2 a;o 

+— A(x ,h(x )) 
x 

= [h'(x )(-2X 2 h(x ) + A 2 + [i + Ai(l - x Q )) - Ai/i(x )]/i 2 a;o 

since A(xq, h(xo)) = 0. Therefore, dA(x, h(x))/dx = at point x = xq iff (note that xo 7^ 0) 

h'{x )(-2X 2 h(x ) + X 2 +fj, + Xi(l-x )) - Xih(x ) = 0. 

Since — 2A 2 /i(xo) + A 2 + \i + Ai(l — xq) < because xo > 1, we get 

Xih(x ) 



h'(x ) 



-2X 2 h(x ) + A 2 + fx + Ai (1 - x ) 



with (see ([78]) /i(x ) = (A + ^i)x /(Ax + ^i), so that h'(x ) < 0. However, h'(x ) > in (J5D]). This yields 
a contradiction, thereby implying that dA(x,h(x))/dx does not vanish at point x = xq when A(x,h(x)) 
does or, equivalently, that xq is a zero of multiplicity one. I 

Lemma A. 4 Under conditions 132) and Assumption A the index \ °f the Riemann-Hilbert problem (the 
index is defined in 1161]) ) is equal to zero. 

Proof. 



Recall the definition of U(x) in ((57)) . First, by studying U(y jli/Xie 1,0 ) for 9 € [0,27r) it is easily seen 
that U(x) describes a closed (and simple) contour when x describes the circle C ^/ -^ I moreover, for 



x G C -i U{x) takes only real values when x € {— y /ti/Ai, y /*i/Ai}. 
As a result, we will show that x = if we show that 



U -VA1/A1 )xU VA1/A1 >0, (82) 
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since ([8"2")l will imply that the contour denned by {U{x) : \x\ = y fii/Xi} does not contain the point x — 
in its interior, so that by definition of the index, x = 0. 

We have from ([5? ]) - ([55 ]) (Hint: R(x,h(x)) = by definition of h(x))) 

A(x,h(x)) = - h(x))x + m{x - h(x)) (83) 

a 

B(x,h(x)) = -^-(X(l-x)h(x) + fj, 2 (h(x)-x)). (84) 
a 



Define X- := — y/ti/Ai and a; + := y /ti/Ai. 

By Assumption A we know that x_ < —1 and x + > 1. Also note that h(x—) = y\ < 1 and /i(.t+) = j/2 < 1 
from Proposition IOI (d2) and (gSJ). With this, it it is easily seen from ([8"3" )) -([8"3 )l that 

A(x-,h(x-)) > and < 

and 

B(x-,h(x-)) < and B(x+,h(x+)) > 

so that 

A(x-,h(x-))/B(x-,h(x-)) < and {A(x+,h(x + ))/B(x + ,h(x+)) < 0. 

and, therefore, 

A(x-,h(x-))/B(x-,h(x-))A(x + ,h(x+))/B(x + ,h(x+)) > 0. (85) 
The above shows that (|52^) is true if r = in the definition ofU(x) since in this case U (x) = A(x, h(x))/B(x, h(x)). 



Assume that r = 1 in the definition of U(x) with x$ < x + and (A + ^i)xo/(Aa;o + A*i) — y A2/A2. Since 
(a; — xo) < for x = x_ and (a; — xq) > for x = x + we conclude from ([85)1 that E/(a;_) > and 
C/(a;+) > 0, thereby showing that (|82p is also true in this case. 



It remains to investigate the case when r = 1 with xq = x+ and (A + /ii)a:o/(A.To + /ii) < y /t2/A2- Clearly, 
U(x-) > since, from ([55)1 . A(x_, h(x-))/B(x-, h(x-)) < and (x_ — Xo) < because £_ < — 1. 

Let us focus on the sign of U{x+). We know that the mapping x — > U(x) is continuous for |x| < 
x + and that C/(x+) ^ when x+ = icq- Since we have shown that U(x+) > when xq < x + and 



(A + H\)X{) I ' (\xq + Hi) < y //2/A2, we deduce, by continuity, that necessarily ?7(x+) > when a* + = 



and (A + /j,i)xq/ (Axo + Mi) < y A2/A2, which concludes the proof. I 

Lemma A. 5 Under condition iSty) and Assumption A, B{k{y) 1 y) = /or |y| = 1, y =/= 1. Afeo, 
B(k(y),y) has a zero at y = 1, uref/i multiplicity one. 

Proof. Fix 1 2/ 1 = 1, 2/ 7^ 1- We know from Proposition 15. 11 - (al) that \k(y)\ < 1. 

From ([3"8"[) and the fact that R(k(y),y) = by definition of k(y), we see that B(k(y), y) = is equivalent 
to 

= A(l - k(y))y + M2 (y - k{y)) = (A(l - fc(y)) + fx 2 )y - fi 2 k(y) 
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that is, 

A(l - k(y) + fj, 2 )y = ^%k(y). 
Taking the absolute value in both sides of the above equation yields 

|A(1 - k(y) +A*a)| = [A(l - k(y) + fi 2 )y\ = \fi 2 k(y)\ < » 2 . (86) 

But |A(1 - k(y)) + /Lt 2 )| > which contradicts (J86J) . Hence, B(k(y),y) ^ for \y\ = 1, y ^ 1. 

Since fc(l) = 1, we see that B(k(l), 1) = B(l, 1) = from the definition of B(x, y). Let us show that the 
multiplicity of this zero is one. This amounts to showing that dB(k(y), y)/dy does not vanish at y = 1. 

Differentiating B(k(y),y) w.r.t. y in (f3"5|) (Hint: R(k(y),y) = 0) and setting y = 1, gives 

d -^^\y=l = ^((X + , 2 )k'(l)-^. (87) 

Let us calculate fc'(l), the derivative of k(y) at y = 1. To this end, let us use (|34[) to differentiate 
R(k(y),y) (which is equal to zero) w.r.t. y, which gives 

= dR ^y) \ = (m _ aAl )fc'(l) + M( u 2 - a\ 2 (88) 

ay 

so that fc'(l) = (aA 2 — MM2) / (Wi — a -^i) (note that — aAi 7^ from Assumption A, which shows 
that fc'(l) is well defined). Plugging this value of k'(l) into (j5Tj) gives 

dB{k{y),y) W*i ,, , \t\ , . \ 1 \ w 

- A \y=i = — f n aA2 A + ^2 -M2 Wi -aAi)) 

dy a(MMi — a Ai) 

= —7 — a(AA2 + A^2 — MM2) 

— aAi) 

MMi / AA 2 A \ 

= + - - 1 < 

^1 - aAi VM^2 M / 

under the conditions in (|32[) (to establish the 2nd equality we have used the definitions of a and A) . This 
proves that dB(k(y), y)/dy\ y= \ 7^ and completes the proof. ■ 
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Figure 2: Transition-rate diagram. 
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Figure 3: Probability of an empty system (fx = 4, fx\ = fj,2 = 2). 




Figure 4: The expected orbit queue size, -E[Q2] (a* = 4, fii = [ii = 2). 
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Figure 5: The expected queue lengths of the orbit queues as functions of fj,2 (Ai = X-2 = 1.2, /x = 4, 
Mi =2). 











EP,] 






E[Q 2 ] 








i 






I 






* 






* 






i* ~ 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 

^2 

Figure 6: The expected queue lengths of the orbit queues as functions of A2 (Ai = 1, fj,i = \i-z = 2, fj, = 4). 
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